Original model:
     f(x,y) = (((((1.3806503E-23.*y)./6.626068E-34).*exp(-(DeltaGCat./(8.314472.*y)
                    ))).*0.000000586)./(((( 1.3806503E-23.*y)./6.626068E-
                    34).*exp(-(DeltaGInact./( 8.314472.*y)))).*(exp((DeltaHEq./8.314472)
                    .*(1./TEq-1./y))))).*(1-exp(-(((((1.3806503E-23.*y)
                    ./6.626068E-34).*exp(-(DeltaGInact./(8.314472.*y)))
                    ).*(exp((DeltaHEq./8.314472).*(1./TEq-1./y))).*x)./(1+
                    (exp((DeltaHEq./8.314472).*(1./TEq-1./y)))))))

Round 1:
Coefficients (with 95% confidence bounds):
       DeltaGCat =  7.819e+004  (7.819e+004, 7.819e+004)
       DeltaGInact =  1.046e+005  (1.046e+005, 1.046e+005)
       DeltaHEq =  1.169e+005  (1.168e+005, 1.17e+005)
       TEq =       346.9  (346.9, 346.9)

Goodness of fit:
  SSE: 1.11e-007
  R-square: 0.9998
  Adjusted R-square: 0.9998
  RMSE: 2.936e-006

Round 2:
Coefficients (with 95% confidence bounds):
       DeltaGCat =  7.819e+004  (7.819e+004, 7.819e+004)
       DeltaGInact =  1.046e+005  (1.046e+005, 1.046e+005)
       DeltaHEq =  1.169e+005  (1.168e+005, 1.17e+005)
       TEq =       346.9  (346.9, 346.9)

Goodness of fit:
  SSE: 1.161e-007
  R-square: 0.9998
  Adjusted R-square: 0.9998
  RMSE: 3.002e-006

Model O:
     f(x,y) = (((((1.3806503E-23.*y)./6.626068E-34).*exp(-(DeltaGCat./(8.314472.*y)
                    ))).*0.000000586)./(((1.3806503E-23.*y)./6.626068E-
                    34).*exp(-(DeltaGDeact./(8.314472.*y))))).*(1-exp(-
                    (((((1.3806503E-23.*y)./6.626068E-34).*exp(-(DeltaGDeact./(8.314472.*y)
                    ))).*x)./(1+exp((DeltaHEq./8.314472).*(1./TEq-1./y)
                    )))))

Round 1:
Coefficients (with 95% confidence bounds):
       DeltaGCat =   7.84e+004  (7.839e+004, 7.84e+004)
       DeltaGDeact =  1.019e+005  (1.019e+005, 1.019e+005)
       DeltaHEq =  1.733e+005  (1.732e+005, 1.735e+005)
       TEq =       350.1  (350.1, 350.1)

Goodness of fit:
  SSE: 2.117e-007
  R-square: 0.9997
  Adjusted R-square: 0.9997
  RMSE: 4.055e-006

Round 2:
Coefficients (with 95% confidence bounds):
       DeltaGCat =   7.84e+004  (7.839e+004, 7.84e+004)
       DeltaGDeact =  1.019e+005  (1.019e+005, 1.019e+005)
       DeltaHEq =  1.733e+005  (1.732e+005, 1.735e+005)
       TEq =       350.1  (350.1, 350.1)

Goodness of fit:
  SSE: 2.116e-007
  R-square: 0.9997
  Adjusted R-square: 0.9997
  RMSE: 4.054e-006

Model P:
     f(x,y) = (((((1.3806503E-23.*y)./6.626068E-34).*exp(-(DeltaGCat./(8.314472.*y)
                    ))).*0.000000586)./((((1.3806503E-23.*y)./6.626068E-
                    34).*exp(-(DeltaGDeact./(8.314472.*y))))+(((1.3806503E-
                    23.*y)./6.626068E-34).*exp(-(DeltaGInact./(8.314472.*y)
                    ))).*exp((DeltaHEq./8.314472).*(1./TEq-1./y)))).*(1-
                    exp(-(((((1.3806503E-23.*y)./6.626068E-34).*exp(-(DeltaGDeact./(8.314472.*y)
                    )))+(((1.3806503E-23.*y)./6.626068E-34).*exp(-(DeltaGInact./(8.314472.*y)
                    ))).*exp((DeltaHEq./8.314472).*(1./TEq-1./y))).*x)./(1+
                    exp((DeltaHEq./8.314472).*(1./TEq-1./y)))))

Round 1:
Coefficients (with 95% confidence bounds):
       DeltaGCat =  7.821e+004  (7.821e+004, 7.821e+004)
       DeltaGDeact =  1.488e+005  (-1.783e+006, 2.08e+006)
       DeltaGInact =  1.047e+005  (1.047e+005, 1.047e+005)
       DeltaHEq =  1.184e+005  (1.182e+005, 1.185e+005)
       TEq =         347  (347, 347)

Goodness of fit:
  SSE: 2.575e-007
  R-square: 0.9996
  Adjusted R-square: 0.9996
  RMSE: 4.472e-006

Round 2:
Coefficients (with 95% confidence bounds):
       DeltaGCat =   7.82e+004  (7.82e+004, 7.82e+004)
       DeltaGDeact =  1.565e+005  (-6.596e+006, 6.908e+006)
       DeltaGInact =  1.046e+005  (1.046e+005, 1.046e+005)
       DeltaHEq =  1.176e+005  (1.175e+005, 1.177e+005)
       TEq =         347  (346.9, 347)

Goodness of fit:
  SSE: 1.663e-007
  R-square: 0.9998
  Adjusted R-square: 0.9998
  RMSE: 3.594e-006

Round 3:
Coefficients (with 95% confidence bounds):
       DeltaGCat =   7.82e+004  (7.819e+004, 7.82e+004)
       DeltaGDeact =   1.53e+005  (-2.11e+006, 2.416e+006)
       DeltaGInact =  1.046e+005  (1.046e+005, 1.046e+005)
       DeltaHEq =  1.172e+005  (1.171e+005, 1.173e+005)
       TEq =       346.9  (346.9, 346.9)

Goodness of fit:
  SSE: 1.131e-007
  R-square: 0.9998
  Adjusted R-square: 0.9998
  RMSE: 2.963e-006

Round 4:
Coefficients (with 95% confidence bounds):
       DeltaGCat =  7.819e+004  (7.819e+004, 7.82e+004)
       DeltaGDeact =  1.567e+005  (-4.679e+006, 4.993e+006)
       DeltaGInact =  1.046e+005  (1.046e+005, 1.046e+005)
       DeltaHEq =  1.171e+005  (1.17e+005, 1.172e+005)
       TEq =       346.9  (346.9, 346.9)

Goodness of fit:
  SSE: 1.214e-007
  R-square: 0.9998
  Adjusted R-square: 0.9998
  RMSE: 3.07e-006

Round 5:
Coefficients (with 95% confidence bounds):
       DeltaGCat =  7.819e+004  (7.819e+004, 7.82e+004)
       DeltaGDeact =  1.568e+005  (-3.923e+006, 4.236e+006)
       DeltaGInact =  1.046e+005  (1.046e+005, 1.046e+005)
       DeltaHEq =  1.172e+005  (1.171e+005, 1.173e+005)
       TEq =       346.9  (346.9, 346.9)

Goodness of fit:
  SSE: 1.171e-007
  R-square: 0.9998
  Adjusted R-square: 0.9998
  RMSE: 3.016e-006

Round 6:
Coefficients (with 95% confidence bounds):
       DeltaGCat =   7.82e+004  (7.819e+004, 7.82e+004)
       DeltaGDeact =  1.592e+005  (-7.842e+006, 8.16e+006)
       DeltaGInact =  1.046e+005  (1.046e+005, 1.046e+005)
       DeltaHEq =  1.173e+005  (1.172e+005, 1.173e+005)
       TEq =       346.9  (346.9, 346.9)

Goodness of fit:
  SSE: 9.917e-008
  R-square: 0.9999
  Adjusted R-square: 0.9999
  RMSE: 2.775e-006

Model M:
     f(x,y) = (((((1.3806503E-23.*y)./6.626068E-34).*exp(-(DeltaGCat./(8.314472.*y)
                    ))).*0.000000586)./((((1.3806503E-23.*y)./6.626068E-
                    34).*exp(-(DeltaGInact./(8.314472.*y)))).*(1+exp((DeltaHEq./8.314472)
                    .*(1./TEq-1./y))))).*(1-exp(-((((1.3806503E-23.*y)./6.626068E-
                    34).*exp(-(DeltaGInact./(8.314472.*y)))).*x)))

Round 1:
Coefficients (with 95% confidence bounds):
       DeltaGCat =  7.828e+004  (7.828e+004, 7.828e+004)
       DeltaGInact =  1.054e+005  (1.054e+005, 1.054e+005)
       DeltaHEq =  1.277e+005  (1.276e+005, 1.277e+005)
       TEq =         348  (348, 348.1)

Goodness of fit:
  SSE: 9.152e-008
  R-square: 0.9999
  Adjusted R-square: 0.9999
  RMSE: 2.666e-006

Round 2:
Coefficients (with 95% confidence bounds):
       DeltaGCat =  7.828e+004  (7.828e+004, 7.828e+004)
       DeltaGInact =  1.054e+005  (1.054e+005, 1.054e+005)
       DeltaHEq =  1.277e+005  (1.276e+005, 1.277e+005)
       TEq =         348  (348, 348.1)

Goodness of fit:
  SSE: 8.273e-008
  R-square: 0.9999
  Adjusted R-square: 0.9999
  RMSE: 2.535e-006

Model N:
     f(x,y) = (((((1.3806503E-23.*y)./6.626068E-34).*exp(-(DeltaGCat./(8.314472.*y)
                    ))).*0.000000586)./(1+(exp((DeltaHEq./8.314472).*(1./TEq-
                    1./y))))).*(x-(((((1.3806503E-23.*y)./6.626068E-34)
                    .*exp(-(DeltaGInact./(8.314472.*y)))).*(((((1.3806503E-
                    23.*y)./6.626068E-34).*exp(-(DeltaGInact./(8.314472.*y)
                    )))+(((1.3806503E-23.*y)./6.626068E-34).*exp(-(DeltaGreact./(8.314472.*y)
                    )))).*x+exp(-((((1.3806503E-23.*y)./6.626068E-34).*exp(-
                    (DeltaGInact./(8.314472.*y))))+(((1.3806503E-23.*y)
                    ./6.626068E-34).*exp(-(DeltaGreact./(8.314472.*y)))
                    )).*x)))./(((((1.3806503E-23.*y)./6.626068E-34).*exp(-
                    (DeltaGInact./(8.314472.*y))))+(((1.3806503E-23.*y)
                    ./6.626068E-34).*exp(-(DeltaGreact./(8.314472.*y)))
                    )).*((((1.3806503E-23.*y)./6.626068E-34).*exp(-(DeltaGInact./(8.314472.*y)
                    )))+(((1.3806503E-23.*y)./6.626068E-34).*exp(-(DeltaGreact./(8.314472.*y)
                    )))))))

Round 1 (robust OFF):
Coefficients (with 95% confidence bounds):
       DeltaGCat =  7.836e+004  (7.795e+004, 7.878e+004)
       DeltaGInact =  9.521e+004  (8.932e+004, 1.011e+005)
       DeltaGreact =  8.945e+004  (8.707e+004, 9.183e+004)
       DeltaHEq =   1.64e+005  (1.631e+005, 1.649e+005)
       TEq =       348.8  (348.6, 349)

Goodness of fit:
  SSE: 2.446e-005
  R-square: 0.966
  Adjusted R-square: 0.966
  RMSE: 4.359e-005

Round 2 (r O):
Coefficients (with 95% confidence bounds):
       DeltaGCat =  7.831e+004  (7.787e+004, 7.876e+004)
       DeltaGInact =  9.469e+004  (8.934e+004, 1e+005)
       DeltaGreact =  8.939e+004  (8.727e+004, 9.151e+004)
       DeltaHEq =  1.651e+005  (1.642e+005, 1.66e+005)
       TEq =       348.9  (348.7, 349.1)

Goodness of fit:
  SSE: 2.413e-005
  R-square: 0.9665
  Adjusted R-square: 0.9665
  RMSE: 4.329e-005

Model Q:
     f(x,y) = (((((1.3806503E-23.*y)./6.626068E-34).*exp(-(DeltaGCat./(8.314472.*y)
                    ))).*0.000000586)./(1+(exp((DeltaHEq./8.314472).*(1./TEq-
                    1./y))))).*(x-(((((((1.3806503E-23.*y)./6.626068E-34)
                    .*exp(-(DeltaGa./(8.314472.*y))))+(((1.3806503E-23.*y)
                    ./6.626068E-34).*exp(-(DeltaGi./(8.314472.*y)))).*(exp((DeltaHEq./8.314472)
                    .*(1./TEq-1./y))))./(1+(exp((DeltaHEq./8.314472).*(1./TEq-
                    1./y))))).*(((((((1.3806503E-23.*y)./6.626068E-34).*exp(-
                    (DeltaGa./(8.314472.*y))))+(((1.3806503E-23.*y)./6.626068E-
                    34).*exp(-(DeltaGi./(8.314472.*y)))).*(exp((DeltaHEq./8.314472)
                    .*(1./TEq-1./y))))./(1+(exp((DeltaHEq./8.314472).*(1./TEq-
                    1./y)))))+(((1.3806503E-23.*y)./6.626068E-34).*exp(-
                    (DeltaGb./(8.314472.*y))))+(((1.3806503E-23.*y)./6.626068E-
                    34).*exp(-(DeltaGj./(8.314472.*y))))).*x+exp(-((((((1.3806503E-
                    23.*y)./6.626068E-34).*exp(-(DeltaGa./(8.314472.*y)
                    )))+(((1.3806503E-23.*y)./6.626068E-34).*exp(-(DeltaGi./(8.314472.*y)
                    ))).*(exp((DeltaHEq./8.314472).*(1./TEq-1./y))))./(1+
                    (exp((DeltaHEq./8.314472).*(1./TEq-1./y)))))+(((1.3806503E-
                    23.*y)./6.626068E-34).*exp(-(DeltaGb./(8.314472.*y)
                    )))+(((1.3806503E-23.*y)./6.626068E-34).*exp(-(DeltaGj./(8.314472.*y)
                    )))).*x)))./(((((((1.3806503E-23.*y)./6.626068E-34)
                    .*exp(-(DeltaGa./(8.314472.*y))))+(((1.3806503E-23.*y)
                    ./6.626068E-34).*exp(-(DeltaGi./(8.314472.*y)))).*(exp((DeltaHEq./8.314472)
                    .*(1./TEq-1./y))))./(1+(exp((DeltaHEq./8.314472).*(1./TEq-
                    1./y)))))+(((1.3806503E-23.*y)./6.626068E-34).*exp(-
                    (DeltaGb./(8.314472.*y))))+(((1.3806503E-23.*y)./6.626068E-
                    34).*exp(-(DeltaGj./(8.314472.*y))))).*((((((1.3806503E-
                    23.*y)./6.626068E-34).*exp(-(DeltaGa./(8.314472.*y)
                    )))+(((1.3806503E-23.*y)./6.626068E-34).*exp(-(DeltaGi./(8.314472.*y)
                    ))).*(exp((DeltaHEq./8.314472).*(1./TEq-1./y))))./(1+
                    (exp((DeltaHEq./8.314472).*(1./TEq-1./y)))))+(((1.3806503E-
                    23.*y)./6.626068E-34).*exp(-(DeltaGb./(8.314472.*y)
                    )))+(((1.3806503E-23.*y)./6.626068E-34).*exp(-(DeltaGj./(8.314472.*y)
                    )))))))

Round 1 (robust OFF):
Coefficients (with 95% confidence bounds):
       DeltaGCat =  7.844e+004  (7.804e+004, 7.885e+004)
       DeltaGa =  9.656e+004  (8.77e+004, 1.054e+005)
       DeltaGb =  9.149e+004  (-1.727e+010, 1.727e+010)
       DeltaGi =  9.367e+004  (4.133e+004, 1.46e+005)
       DeltaGj =  9.149e+004  (-1.727e+010, 1.727e+010)
       DeltaHEq =  1.678e+005  (1.636e+005, 1.721e+005)
       TEq =       349.6  (329.3, 369.9)

Goodness of fit:
  SSE: 2.364e-005
  R-square: 0.9672
  Adjusted R-square: 0.9672
  RMSE: 4.286e-005

Round 2 (r O):
Coefficients (with 95% confidence bounds):
       DeltaGCat =  7.825e+004  (7.753e+004, 7.898e+004)
       DeltaGa =  9.322e+004  (8.615e+004, 1.003e+005)
       DeltaGb =  9.048e+004  (-1.433e+010, 1.433e+010)
       DeltaGi =  9.319e+004  (-8212, 1.946e+005)
       DeltaGj =  9.048e+004  (-1.433e+010, 1.433e+010)
       DeltaHEq =  1.696e+005  (1.654e+005, 1.738e+005)
       TEq =       349.2  (315.2, 383.2)

Goodness of fit:
  SSE: 2.331e-005
  R-square: 0.9676
  Adjusted R-square: 0.9676
  RMSE: 4.255e-005

Round 1 (lar on): compared with Round 1 (robust OFF) above:
Coefficients (with 95% confidence bounds):
       DeltaGCat =  7.676e+004  (7.317e+004, 8.035e+004) less certain, close tho
       DeltaGa =  8.447e+004  (7.585e+004, 9.31e+004) ok
       DeltaGb =  8.423e+004  (8.238e+004, 8.608e+004) much more certain
       DeltaGi =  8.729e+004  (-3.33e+004, 2.079e+005) less certain
       DeltaGj =  1.329e+005  (-2.574e+009, 2.574e+009) less certain and rubbish
       DeltaHEq =  1.594e+005  (1.504e+005, 1.683e+005) close
       TEq =       346.4  (276.9, 415.9) a bit off, lower, less certain
Goodness of fit:
  SSE: 1.899e-006
  R-square: 0.9974
  Adjusted R-square: 0.9974
  RMSE: 1.215e-005

Round 2 (lar on):
Coefficients (with 95% confidence bounds):
       DeltaGCat =  7.683e+004  (7.239e+004, 8.128e+004)
       DeltaGa =  8.447e+004  (7.334e+004, 9.56e+004)
       DeltaGb =   8.41e+004  (8.146e+004, 8.673e+004)
       DeltaGi =  8.733e+004  (-6.799e+004, 2.426e+005)
       DeltaGj =  1.354e+005  (-7.349e+009, 7.349e+009)
       DeltaHEq =  1.606e+005  (1.52e+005, 1.691e+005)
       TEq =       346.5  (262.7, 430.4)

Goodness of fit:
  SSE: 3.213e-006
  R-square: 0.9955
  Adjusted R-square: 0.9955
  RMSE: 1.58e-005

Round 1 (robust off) parameters into a lar fit:
didn't converge. Lar doesn't give sensible parameters yet gives stats, hence I don't trust it with a barge-pole.

Model R:
     f(x,y) = (((((1.3806503E-23.*y)./6.626068E-34).*exp(-(DeltaGCat./(8.314472.*y)
                    ))).*0.000000586)./(1+(exp((DeltaHEq./8.314472).*(1./TEq-
                    1./y))))).*(x-((((((1.3806503E-23.*y)./6.626068E-34)
                    .*exp(-(DeltaGInact./(8.314472.*y))))./(1+1./(exp((DeltaHEq./8.314472)
                    .*(1./TEq-1./y))))).*((((((1.3806503E-23.*y)./6.626068E-
                    34).*exp(-(DeltaGInact./(8.314472.*y))))./(1+1./(exp((DeltaHEq./8.314472)
                    .*(1./TEq-1./y)))))+(((1.3806503E-23.*y)./6.626068E-
                    34).*exp(-(DeltaGreact./(8.314472.*y))))).*x+exp(-(((((1.3806503E-
                    23.*y)./6.626068E-34).*exp(-(DeltaGInact./(8.314472.*y)
                    )))./(1+1./(exp((DeltaHEq./8.314472).*(1./TEq-1./y)
                    ))))+(((1.3806503E-23.*y)./6.626068E-34).*exp(-(DeltaGreact./(8.314472.*y)
                    )))).*x)))./((((((1.3806503E-23.*y)./6.626068E-34).*exp(-
                    (DeltaGInact./(8.314472.*y))))./(1+1./(exp((DeltaHEq./8.314472)
                    .*(1./TEq-1./y)))))+(((1.3806503E-23.*y)./6.626068E-
                    34).*exp(-(DeltaGreact./(8.314472.*y))))).*(((((1.3806503E-
                    23.*y)./6.626068E-34).*exp(-(DeltaGInact./(8.314472.*y)
                    )))./(1+1./(exp((DeltaHEq./8.314472).*(1./TEq-1./y)
                    ))))+(((1.3806503E-23.*y)./6.626068E-34).*exp(-(DeltaGreact./(8.314472.*y)
                    )))))))

Round 1 (lar):
Coefficients (with 95% confidence bounds):
       DeltaGCat =  7.846e+004  (7.846e+004, 7.847e+004)
       DeltaGInact =  1.102e+005  (1.079e+005, 1.125e+005)
       DeltaGreact =  1.002e+005  (9.917e+004, 1.012e+005)
       DeltaHEq =   1.54e+005  (1.537e+005, 1.544e+005)
       TEq =       347.9  (347.8, 347.9)

Goodness of fit:
  SSE: 1.276e-006
  R-square: 0.9982
  Adjusted R-square: 0.9982
  RMSE: 9.955e-006 ----- didn't converge with 1000 FUNC EVALS but sensible as parameters!

Round 2 (lar) with 10000 func evals, did convergre:
Coefficients (with 95% confidence bounds):
       DeltaGCat =  7.846e+004  (7.845e+004, 7.846e+004)
       DeltaGInact =  1.102e+005  (1.079e+005, 1.124e+005)
       DeltaGreact =  1.002e+005  (9.915e+004, 1.012e+005)
       DeltaHEq =  1.537e+005  (1.534e+005, 1.54e+005)
       TEq =       347.8  (347.7, 347.9)

Goodness of fit:
  SSE: 1.257e-006
  R-square: 0.9983
  Adjusted R-square: 0.9983
  RMSE: 9.879e-006